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Abstract 



q} . We study the behaviour of two Brownian particles coupled by an elastic 

harmonic force in a quenched disordered medium. We found that to first 
order in disorder strength, the relative motion weakens (with respect to the 
reference state of a Brownian particle with the double mass) the effect of 
the quenched forces on the centre of mass motion of the Brownian particles, 
so that the motion will become less subdiffusive (superdiffusive) for potential 
(solenoidal) disorder. The mean-square relative distance between the particles 
t^ \ behaves in a different way depending of whether the particles are free to move 

or one particle is anchored in the space. While the effect of nonpotential 
disorder consists in increasing the mean-square distance in both cases, the 
potential disorder decreases the mean-square distance, when the particles are 

free to move, and increases it when one particle is anchored in the space. 
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Diffusion of a particle in quenched random environment has been the subject of numerous 
studies in the past decade [|J- [|10| . It is of interest in problems like tracer diffusion in porous 
media, transport in nonsymmetrical hopping systems, diffusion in a fluid with stationary 
random streams. For short range disorder the diffusion is anomalous, < x 2 (t) > ~ t 2u with 
v 7^ 1/2, for space dimensions lower than 2. In this Letter we present the results of the 
study of the behaviour of two Brownian particles inserted in the same disordered media and 
coupled by a harmonic force. The main interest in this problem is due to the problem of 
the behaviour of a polymer chain in disordered media [|ll|]- flnj. Two harmonically coupled 



Brownian particles in disordered media may be considered as the simplest case of a polymer 
chain, which has only two modes: the centre of mass mode and the relative mode, and 
thus may serve as a toy model of the much more complicated problem of a polymer in 
disordered media. We will study in the present Letter the interplay between the centre of 
mass motion and the relative motion in the presence of quenched disorder. Recently, the 
interplay between the centre of mass motion and the relative motion has been considered 
in the case of two electrons in disordered media [|19|]- | 22| . The task we study here is the 



classical counterpart of the problem of two coupled electrons in disordered media |19|- |22 
The Langevin equations for the coupled Brownian particles are 

r y-Tjr + X(xi-x 2 )=ri 1 {t) + F{xi), 1~£ + A(x 2 - x 1 ) =r? 2 (*) + F(x 2 ), (1) 

where 7 is the monomer friction coefficient, A is the spring constant, 7/1 (£) is the thermal noise, 

which is assumed to be Gaussian distributed with zero mean and the correlator Vi(t)r]{(t') = 
2 / ykT5 l: >5(t — t') and fcT/7 = Do being the diffusion constant of the Brownian particle. The 
quenched random force F(x) is assumed to be Gaussian distributed with the zero mean and 
the correlator < F l (x)F : >(x / ) > = C^(x — x'), where the Fourier transform of C^(x — x') is 
given by 

C^(q) = q-^iAqttf/q 2 + B{5 lj - q l q j /q 2 )). (2) 



By using the Martin-Siggia formalism [2^] the Langevin equations (|1|) can be represented 
as a path integral over the trajectories of the Brownian particles in the phase space, which 
is composed of the coordinates of the particles, x a (t), and the response fields (momenta), 
Pa{t), with the weight determined by the action S = So + Si with So and Si given by 

S = kT-f J dt'p 2 a (t')-i f dt'p a {t'){^dx a {t')/dt' -iX I ' dt'p l {t'){x l {t')-x 2 {t')) - 
t\Jdt'p 2 (t')(x 2 (t')-x l (t'))), (3) 

Si = - dti dt' 2 p a (ti)C a/ 3(x a (ti) -xp(t 2 ))pp(t 2 ), (4) 

where the sum convention over Greek indices is used above. 

To study the perturbation expansions of the quantities such as the centre of mass and 
the relative distance between the particles we consider the correlation function 



< exp(p c x c (t) + ip r x r (t)) > , (5) 

where the coordinate of the centre of mass x c and the relative coordinate x r are defined as 
follows x c = (x\ + X2V2, and x r = X\ — X2- This transformation of the coordinates implies 
the transformation of the response fields in ([3]) and momenta in (|5D as follows p c = Pi + Pi 
and p r = (pi — p 2 )/2. The subscript in Eq.(|^) means that x c (t ) and x r (t ) are fixed at the 
initial time to- The average over x r (to) with the equilibrium Boltzmann weight is obtained 
by taking the limit t — ► —00. In studying < x c {t) 2 > we consider x c (to) to be fixed (= 0), 
while in computing < x r (t) 2 > we proceed as follows. We confine the centre of mass mode 
in a harmonic potential with an elastic constant Ao, let to tend to —00 and then put Ao = 0. 
This procedure corresponds to the integration over x c (t ) in volume ~ (D /\ ) d ^ 2 . Due to 
the procedures of carrying out the averages over the initial conditions the results we obtain 
below, Eqs.(|| [10]), are restricted from below by a small nonzero value of the harmonic 
coupling A. 

The average in Eq.(|5]) can be represented as path integral over the phase trajectories of 
Brownian particles weighted with MSR functional (|3|-(|). Due to introduction of the centre 
of mass and relative coordinates the free part of the action S splits respectively into two 
independent parts dependent of the centre of mass and the relative coordinates, which, 
however, are coupled due to the interaction part of the action Si. 

The mean-square displacement < x c (t) 2 > of the centre of mass is obtained by using ([5]) 
to first order in powers of the disorder strength as 



< x c (t) 2 > = 2dD c t + ^|p!W2 - x)D; d ^J^ 



dy(t - y) x 



{(A + dB- B)(T(y)- d/2+ " + U(y)- d/2+ *) - 2Ay(d/2 - x) x 

{ny) -i- d/ 2 + « + u{y) -i~ d /2 + « _ eM -\y)(U{y)- l - d / 2+ " - Tiy)- 1 -^}, (6) 

where T(y) = y + (1 - exp(-Ay))/(A), U{y) = y + (1 + exp(-Ay))/(A) with A = 2A/ 7 , 
D c = -Do/2- In the limit of large A we obtain from Eq.(^) the first-order correction to 
the mean-square displacement of a Brownian particle with the diffusion constant D c . After 
carrying out the integration over time we obtain 

1 9 
< x {tf > = 2dD c t(l + ——-L-T(d/2 - x)D~ 2 2(A + dB - B - 2A(d/2 - x)) x 

2d2(2ii) d 

(D c t) l+ *- d / 2 + ...). (7) 



(1 + x-d/2){2 + x-d/2) 

We have found that the preasymptotic term in the expansion of Eq. (^|) for large A does not 
possess the 1/e-pole (e = d c — d) in the vicinity of the upper critical dimension d c = 2 + 2k. 
Due to the nonexistence of the fixed-point for the effective coupling constant for potential 
disorder quantitative results beyond the first-order perturbation theory are hard to be made. 
In Fig.l we present the results of the numerical evaluation of the integral in Eq.(^) for 
dimensionalities d — 1.5 and d — 1. For potential field the correction decreases in both 
cases with increasing A. This can be explained qualitatively as follows. When one particle 
is under the action of a well, its partner, which is outside the well will diminish the effect 
of the well, with the result that the first order correction will become less subdiffusive. The 



qualitative run of lines of the force is shown in Fig. 2a. The total effect of the disorder on 
the behavior of the particle is hardly to be predicted, because for potential disorder the 
effective coupling increases under renormalization 0- M. Fig. lb shows the dependence on 
A for d = 1. The first-order correction is positive and decreases to zero with increase of A, 
in agreement with the fact that the first-order correction is zero for one Brownian particle 
@- |§. In this case the partner contribute to average out the effect of the barriers and wells 
acting both subdiffusive, so that for finite A the correction is positive (superdiffusive). A 
prediction of the whole effect of the disorder cannot be made by the present study. 

The reason for superdiffusion in the solenoidal field is due to the remnant of the ballistic 
motion of the particle along the lines of the force (see Fig. 2b), which is disturbed by the 
thermal noise. For two particles the effect is expected to be maximal, if both particles follows 
the same line of the force, which takes place, if the particles are tied closely together (large 
A). Thus, for intermediate A the superdiffusive effect of the solenoidal field will be reduced 
in agreement with the behaviour shown in Fig. la. 

We now will consider the relative motion of two Brownian particles in the case when 
one of the particles is anchored in the space. This problem is equivalent to the case of one 
particle subjected to the external harmonic potential and quenched random forces. The 
mean-square distance between the particles is obtained as 



2 2dD S d Y(d/2-x) 
< x r >= — — h 



2A 2(27r) d A 
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l-exp(-A 2/ ))T 1 (y)- 1 - rf /^}, (8) 



where Ti(y) = (1 — exp(— Xy))/\ and A in (f|j9|) is measured in units of I/7. Carrying out 
the integrations over y and q yields 

Eq.(||) shows that the first-order correction is positive for disorder given by (0), so that the 
effect of disorder consists in increasing the distance between the particles. The increase of 
the relative distance between the particles can be interpreted as the decrease of the elastic 
constant. In the case A = B and k = this result was obtained previously in [23]. The 



increase of the mean-square distance between the particles can be understood qualitatively 
as follows. The interplay between the harmonic potential and the quenched disorder becomes 
important on distances where the harmonic potential is comparable to kT, r 2 ~ D /X. The 
distances with r > y/D /X, which are practically unaccessible in the absence of disorder, 
may become accessible, if the quenched random forces compensate the harmonic potential. 
Due to this the mean-square distance from the origin or the relative distance between the 
particles increases. Notice that according to Eq.([|) the first-order correction for potential 
disorder (B = 0) is regular at the critical dimension, 1 — d c /2 + x = 0, while for solenoidal 
disorder (A = 0) the first-order correction diverges at the critical dimension d c . 

We now will consider the first-order correction to the mean-square distance between the 
Brownian particles in the limit t —>■ 00 by taking into account the diffusion of the centre of 
mass. In computing linu^oo < x 2 (t) > we confine the centre of mass mode in a harmonic 
potential with a small elastic constant Aq, carry out the limit to ~~ *• — °° an d then put 



Ao = 0. This procedure corresponds to the integration of the centre of mass over the volume 
(D /\ ) d / 2 . As a result we have obtained 



€ r 2 ^ _ 2dD o (1 , S d T(d/2- k) ,jsi-d/2+Kr a. 



poo 

2 d/2-«-i {d _ l)B I dyexp(-y)(T- x=1 (y)- d ^ - f/ I=1 (l/)^ /2+ -)} + ...)■ (10) 

Jo 

The integral over y in Eq. (|T0| ) can be tabulated for particular values of d and x. For d = 1 
and x = it is equal to 0.669. It increases with increasing d. It follows from Eq . (|i~0|) that 
the potential disorder (B = 0) results in a decrease of the mean-square relative distance 
between the particles. As it follows from Eq.(|9|) this is in contrast to the case when one 
particle is anchored in the space. 

Let us discuss now the predictions of Eq.([10|) qualitatively. For solenoidal field the 
correction to < x 2 > is positive for both cases with and without the motion of the centre 
of mass. The solenoidal field, which does not possess sources and sinks acts rather as an 
additional thermal noise, which consequently will increase the distance between the particles. 
In the case of the potential disorder the effect with moving centre of mass is opposite to that 
discussed above when one particle is anchored. The decrease of the relative distance can be 
explained to be due to the fact that both particles will prefer to be in the regions where the 



potential has low values. A similar behaviour was predicted by Cates and Ball |13| to occur 
for a polymer chain in disordered media. The decrease of the effect of the wells on the centre 
of mass motion due to the harmonic coupling between the particles is consistent with the 
decrease of the relative distance. Notice that the correction is regular at the upper critical 
dimension, so that at this order the effect is small, and is not expected to give a rise for a 
power law for the effective elastic constant. The decrease of the distance is in agreement 
with the prediction made in [16| for polymer chain, however these authors do not state the 



difference between potential and solenoidal disorder. To our knowledge the present work 
is the first one where the difference in the behaviour with and without the motion of the 
centre of mass is demonstrated in a straightforward calculation. The first-order correction 
to the mean-square relative distance is however regular at the critical dimension for both 
cases associated with Eqs.® and ([K]). The first-order correction to < x 2 > for solenoidal 
disorder (A = 0) is positive and contains a 1/e-pole in the vicinity of the critical dimension. 
To conclude, we have considered the behaviour of two Brownian particles coupled by an 
harmonic potential in a disordered medium to first order of the perturbation theory. We have 
found that the harmonic coupling weakens (with respect to the reference state of a Brownian 
particle with the double mass) the effect of the quenched forces on the centre of mass motion 
of the Brownian particles, i.e. the relative motion contribute in averaging out the effect of 
the disorder. For solenoidal (potential) disorder, the centre of mass motion will become 
less super diffusive (sub diffusive). The mean-square relative distance between the particles 
behaves in a different way depending on whether the Brownian particles are free to move 
or one particle is anchored in the space. While the effect of solenoidal disorder is the same 
in both cases, and consists in increasing the mean-square distance, the potential disorder 
decreases the mean square distance, when the particles are free to move, and increases it 
when one particle is anchored in the space. 
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Figure captions 

Fig.l Dependence of the correction to < x c (t) 2 > as function of A. a) d — 1.5, k = 0; b) 
d — 1, k — 0. 

Fig. 2 The lines of force of a) the potential and b) solenoidal quenched random field. 
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Fig. 1a 
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